Friedel oscillations in quantum degenerate collisional plasma. 

Screening of point charge 



Research of influence of collisions on Friedel oscillations in quantum degen- 
erate collisional plasma (T = 0) is carried out for the first time. It is shown that 
presence of collisions in plasma leads to exponential decreasing of amplitude and 
phase shift of Friedel oscillations. In linear approximation the phase shift is equal 
to the half of quantity inverse to product of Fermi's wave number by free length 
path of electrons. The correct expression for longitudinal dielectric permeabil- 
ity of the quantum collisional plasma found by the authors (see arXiv: 1001.3937 
[math-ph] 22 Jan 2010) is used. 
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Under classical consideration of degenerate electronic plasma poten- 
tial V(r) around point charge Q is described by the classical formula of 
Thomas — Fermi screening [T] 
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1. Introduction 



V(r) = — exp(— kxFr). 
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Here krF is the inverse screening radius of Thomas - - Fermi 

/67me 2 \ V 2 



k TF = 



where e is the electron charge, n is the electron concentration, e F is the 
Fermi energy. 

Fridel was the first [2] [6] who has found out that asymptotic (on 
the large distances) decreasing of screening potential of point charge 
under quantum consideration of degenerate collisionless plasma has not 
only monotonously decreasing, but also oscillatory character. The reason 
of such oscillations is sharp falling (to zero) of Fermi distributions for 
electrons Jf{v) behind Fermi's surface , 

f F (y) = 0{V F - V), 

where Q(x) is the Heaviside function, 



G(x) 



This singularity of Fermi distribution leads to so called Kohn singu- 
larities (see [7j [13] )• Kohn singularities are consequence of logarith- 
mic singularities of the longitudinal dielectric permeability of degenerate 
plasma. Just Kohn singularities lead us to Friedel oscillations in degen- 
erate plasma. 

In the work [TB] the authors found out the dependence of Friedel's 




oscillations on temperature in collisionless non-degenerate plasma. It is 
shown that under non-zero temperature Friedel oscillations amplitude 
decreases exponentially with distance. 

In the present work the research of influence of collisions on Friedel 
oscillations in quantum degenerate collisional plasma (T = 0) is carried 
out. The correct expression for longitudinal dielectric permeability of 
the quantum collisional plasma found by authors (see [H]) is used. It 
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is shown that the presence of collisions in plasma results in exponential 
decreasing of amplitude and phase shift of Friedel oscillations. 

In linear approximation phase shift is equal to half of the quantity 
inverse to product of the Fermi wave number by free length path of elec- 
trons. In more details, phase shift is proportional to effective frequency 
of electron collisions with plasma particles. From here it is clear that 
when the effective frequency of collisions decreases to zero (plasma be- 
comes collisionless), phase shift tends to zero. It is agreed with classical 
result HH, H2J. 

Let's emphasize that we use the expression for dielectric longitudinal 
permeability of collisional degenerate plasma found by us in [14] . This 
expression is deduced on the basis of the solution of the quantum kinetic 
equation for Wigner function in coordinate space, instead of momentum 
space as Mernin [15]. Mermin in [15] has received the expression for 
longitudinal dielectric permeability in the quantum collisional plasma 
using the kinetic equation in momentum space. However, if we use 
the Mermin expression for longitudinal permeability, then exponential 
decreasing of amplitude of Friedel oscillations cannot be found out. As 
the reason of this fact is that the Mermin longitudinal permeability for 
a static case does not depend on collision frequency. 



Let us consider a point charge Q in the origin of coordinates. Then 
the equation describing behaviour of potential V round the point charge 
has the form 



Here A is the Laplace operator, 5(r) is the Dirac delta function of 
vector argument, 5(r) = 5(x)5(y)5(z); p(r) is the density of induced 
charge. 



2. Potential of point charge 




(i.i) 
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We use representation of quantities in (1.1) in the form of Fourier 
integrals 

p(r) = jhf J e ' kI p kd3k - 

Then from the equation (1.1) we receive 

fcV k = 47rp k + AttQ. (1.2) 
Earlier we have shown [141 that 



kj k <7 Z (k)kE k iai{k)k 2 
Pk = , Pk = = V^k- 

CO CO CO 

Here cr^(k) is the longitudinel electric conductivity of quantum plasma 
found in work [H] . 

Substituting this relation into (1.2) we receive 



or 

,2 



^(1 + 1^)^ = 4.0, 



k z e,(k)V k = inQ, 

where ei(\z) is the longitudinal dielectric permeability of quantum plasma, 

sM = i + 4 ™ i(k) 



Thus spectral density of electric potential is expressed in terms of the 
quantity of charge 

When deducing (1.3) we used the fact that dielectric permeability 
depends only on the module k, £/(k) = £i(k). It is obvious that the 
following relation relation V(r) = V(r) also takes place. 
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Then for the quantity V{r) we have 

4ttQ f e ikr 47rQ f e ikrcosd 
V(r) = . xo / — — T-rrd^k = , , n / — — sin 6d6dk. 

Integrating by the angular variable 9 we receive 



X 



N 1 2Q f sin(kr) 

v(r) = - r ■ - J imi dk - (L4) 



o 



The formula (1.4) corresponds to the formula received in the mono- 
graph of Harrison (see | jlll |. p. 334) if we present 



o Q 



q 47TQ 2 ' 

Let's consider asymptotic behaviour of potential V(r), r — > oo. It is 
possible to consider electric field to be weak enough in this area. There- 
fore linear approximation is applicable here. Fridel oscsillations appear 
because of logarithmic singularities of dielectric permeability. 

After double integration of expression (1.4) by parts we find |TT] 

oo 

v(r) = LMf g^MW dfc + (e-*n, 



r 3 7T J kef(k) 



or, believing k = qkp. 



v^y^J^^ dq+ o { ^n. (i.5) 

o 

In (1.5) we have left in explicit form only that term which results 
in oscillations of potential. Other item designated as 0(e~ kTFr ) quickly 
decreases in exponential way, on the Thomas - - Fermi radius it decays 
with distance. 

Integrals of such kind as (1.5) usually tend to zero at tendency of 
distance r to infinity because of fast oscillations sin rk. However as it will 
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be seen the second derivative of e'({k) contains singular Cauchy kernel 
that gives the non-zero contribution to integral. Friedel oscillations are 
caused by logarithmic singularity of dielectric permeability. 

Function standing under the sign of integral in (1.5) is even, therefore 
we can expand the integration onto whole real axis: 

oo 

V[r) = J_ . Q f E^MM rf, + (e-^n, r -» oo. (1.6) 



r° 7r 



qef(q) 



Let us recall that under positive values of q > the quantity q equals 
|k| 

to q = - — , and under negative values the quantity q has not any physical 

ky 

meaning. 

3. Friedel oscillations 

In the expression (1.6) there is an expression of the longitudinal di- 
electric permeability of quantum degenerate collisional plasma £i(q). Ac- 
cording to the work [H] this expression has the following form for statical 
case to = considered here 

ei («)=i+S 1 T (g)+ A (g) - ^ 

2q z 1 - iyg {q) 
In (2.1)the following dimensionless parameters are introduced 



LUr, V 1 



x P = t—^, V = r~ — = 77—, I = vft, 



where u p is the plasma (Langmuir) frequency, hp = Pf/^ is the Fermi 
wave number, pp = mvp is the Fermi momentum, vp is the Fermi veloc- 
ity, m is the electron mass, v is the effective electron frequency, / = vpjv 
is the electron free length path, r = 1/V is the time of electron free length 
path. 

These functions go(q) , g+{q) , g~{q) in statical limit (u — > 0) have the 
following form 

9m) = 7r ln - ' 

2q ty-q 
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g+(q) = ^[(q 2 + ^vf - V 
9-(q) = ^[(q 2 - ^vf - *t 



In 
In 



q 2 + 2q + 2n/ 
q 2 — 2q + 2iy ' 
q 2 -2q- 2iy 
q 2 + 2q-2iy~ 



Let's find roots of the equations 

± 2g ± 2iy = 0. 

We have 

9i,2 = Tl ± V 1 ± 2*j/. (2.2) 

Further let's consider the quantity ?/ as small parameter and designate 
it through £ 

e = y 



v 



kpvp mvp 2e 



In linear approximation by e for roots of (2.2) we receive 

gi,2 = Tl ± (l±ie). 

We present the functions g+(q) and g~(q) in the following form 

{q + ie){q + 2-ie) 
(q — ie)(q — 2 + is)' 

(q + ie)(q — 2 — is) 



2 2 
2 2 

<?-(<?) = ^-(q - 2 - «0(« + 2 + fc) ln 



or, in linear approximation, rejecting terms proportional to e 2 , we have 



9+(q) = ^9i(q); 



g-(q) = ^92(q): 



where 

9l (q) = (q-2 + ie)f+-(q) - (q + 2 - ie)f-+(q) 
g 2 {q) = (q + 2 + - (g - 2 - 

Here the functions containing Kohn singularities are introduced 



f++(q) = (q + 2 + ie) \n(q + 2 + ie), 
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f-(q) = (q-2-ie)Hq-2-i £ ), 
/+-(«) = (g + 2 - ie) ln(g + 2 - ie), 
/-+(?) = (Q-2 + ie)ln(g-2 + ie). 
The second derivatives of these functions result in Cauchy kernels 

Let's find the second derivative e'((q). In the expression for e'({q) we 
leave only those items, which contain Kohn singularities and result in 
the Friedel oscillations 

, = K &ld-J-b) = H 9'!(q) - g'M (9 ^ 

l[Q) 2q 2 l-9o(q) 16g 3 1 - 9 o(q) * 1 ^ 

Let's return to the integral (1.6) and by means of (2.3) we present it 
in the form 

v[r) = i . s M [ m -am y + 0{e - kTni r _ oc . 



r 



16tt J q*e 2 (q)[l-g (q)} 



— oo 



(2.4) 

Let's show that in a statical limit (u — > 0) the expression £i(q) is real. 
Really, according to the results from [H] the dielectric permeability is 
expressed by equality 

f (1 - t 2 )dt 



£i(q) = i + 



3x 



2 J (l-iUT + iqt/y) 2 + q A /(4y 2 ) 



p 



4^ 



1 ) dt 

2 J 1 — iur + zqtf/y 



-l 

We receive from here in the statical limit 

/ (1 - t 2 )dt 



£i(q) = i + 



0^2 J {iqt + y) 2 + q y4 

1 _y I dt 



iqt + 2/ 
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After simple transformations we reduce the previous expression to the 
form which does not contain imaginary unit 

1 



eM) = 1 + 



3i£ p 



{l-t 2 ){qy3 + y 2 -qH 2 )dt 
{q A /A + y 2 -q 2 t 2 ) 2 + Aq 2 y 2 t 2 



i- y - 

2 



/ 

-l 



dt 



q 2 t 2 + y 2 



In the expression g'{ (q) — g'^q) from the equality (2.4) we will leave the 
terms containing Kohn singularities and leading to Friedel oscillations 
again. As the result we receive 



q - 2 + is 



9"{q) - g'lio) = 

q + 2-is q + 2 + is 



+ 



Q 



IE 



q + 2 — ie q — 2 + ie q — 2 — is q + 2 + iE 
Now integral (2.4) we will present in the following form 



1 SQx 



sm(kFfq) 



16tt J q*e*(q)[l-go(q)] 



q-2 + is q + 2 



IE 



— oo 



q + 2 + ie q 
+ 



q + 2 — ie q — 2 + iE 
2 



IE 



q — 2 — ie q + 2 + is 



dq. 



(2.5) 



Let's designate now 



= q+(2 + iE) q + (2 - ie) 
^ W q-(2 + ie)^ q-(2-ie)' 

It is obvious that other two terms in the square brackets are equal to 



<p(-q) = 



q ~ (2 + ie) + q - (2 - ie) 



q+{2 + iE) ' q-\-(2-ie)' 
Now the expression (2.5) we will write in the form 



V( \ - — ^Q x p 
[r) ~ r 3 ' 16tt 



<p(q) - if(-q) 
<?tf(q)[l - 9o(q)] 



sm(kprq) dq. 



— oo 
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This equality can be simplified 



cp(q) sin(kFrq) 
8tt J q^(q)[l-g (q)} 



dq. 



■2.6) 



Let's notice that expression q A £f(q) does not vanish in the point q = 0. 
Besides, we will notice that in the linear approximation of 1— go(q) = 
1. Really, we multiplay the subintegral function from expression for 
go(q) by fraction, numerator and denominator of which is the expression, 
conjugating to the denominator from go(q), we receive: 

l 



9o(q) = ^ 



-i 



(qt + iy) dt 
q 2 t 2 + £ 2 



tyq 



tdt 



dt 



l l 

+ - / 

/ q 2 t 2 + e 2 2 J q 2 t 2 + e 5 

-l -l 



0(e* 



£ 



0), 



as it was required to show. 

Now the integral from (2.6) we will write down in an explicit form 



V(r) 



r 



1 3Qx 

~3 ' 



P 



1 



8 2m 



q+(2 + i£) + q+(2 



l£ 



q 



2 + l£ 



q 



l£ 



9^(9) 



dq. 



-oo 



In last integral subintegral expression has two singular Cauchy ker- 
nels, equal to infinity in the complex conjugate points q = 2±ie. Integra- 
tion near the point q = 2 gives the basic contribution to the expression 
for potential. For the approximate calculus of last integral we use the 
method stated, for example, in the monograph [11]. 

For this purpose we calculate continuous functions from subintegral 
expression in the point q = 2: 



q%(q) 



q=2 



= 16(1 + ^ 2 
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We receive that 



V(r) 



1 3Qxp 
^2(8 + 3x2) 2 



X 



X 



q + (2 + is) q + (2 - iej 



- (2 + ie) g - (2 - ie) 



dq. 



-oo 



The integral from the previous equality is equal to the sum of residues 
relatively the simple poles in points q = 2 ± is: 



3Qx~ 



V ^ r 3 ' (8 + 3x2)2 
From here we obtain 



[2 + ie)e lkFr ^ +l£) + (2 - ie)e 



ikpr(2—ie) 



1 6Qx 2 p e~ kFr£ 
V ^ = ^ ' (8 + 3x?) 2 



2 cos 2/ci?r — e sin 2/c^r 



Finally we get 



-kpre 



V(r) = A — — cos(2k F r + (p), 



or 



V(r) = A 6 -^— cos(2k F r + (2.9) 
r <j 2 



Here 



e e y 
^ = arctg- = - = -, 

12Q^_ 

(8 + 3x2)2- 

The expression (2.9) can be rewritten in dimensional parameters 

v 

e / y 

V(r) = A 5 — cos ( 2k F r + 



r J \ 2kpVF' 

From this formula at v = the classical result for collisionless plasma 
is received (see, for example, fTT] ) 

V(r) = -^-cos(2/ci?r). 
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Let's designate through R = kpr dimensionless length, and introduce 

dimensionless frequency of collisions y = e = . Then the formula 

(2.9) for potential will be written in the form 

o er Ry y 
V(R) = Ak F —3- cos(2i? + |). (2.10) 

Graphical study of Friedel oscillations we will carry out for the case, 
when Akp = 10 5 (fig. 1 - fig. 3 see). From fig. 3 it is visible that 
with growth of frequencies of electron collisions the amplitude of Friedel 
oscillations decreases. 

4. Conclusion 

In the present work we use the expression for dielectric permeability 
of quantum degenerate collisional plasma found by the authors earlier 
[H] . With the help of it the research of influence of collisions on Friedel 
oscillations is carried out for the first time. 

It is shown that presence of collisions in plasma results in exponen- 
tial decrease of amplitude of Friedel oscillations and to inverse shift of 
the phase oscillations. The logarithmic decrement of decrease in linear 

Up 

approximation is equal to r = — , and the inverse phase shift is equal 

v 

v 1 

to — = — — -, where / is the mean free path of electrons. 

2kpVF 2kpl 




Figure 2: Friedel oscillations in the case y = 10 2 , 50 < R < 100. 



14 




15 



References 

1. D. Pines. Elementary exitations in solids. W. A. Benjamin, INC, 
New York - Amsterdam. 1963. 

2. J. Friedel. The distribution of electrons round impurities in monova- 
lent metals. - Phil. Mag. 1952. V. 43, 153 - 189. 

3. J. Friedel. Structure electronique des impuretes dans metaux. - Ann. 
Phys. Paris. 1954. V. 9, 158 - 202. 

4. J. Friedel. Electronic structure of primery solid solutions in metals. 
- Advan. Phys. 1954. V. 3, 446 - 507. 

5. J. Friedel. On some electrical and magnetic properties of metallic 
solid solutions. - Can. J. Phys. 1956. V. 34, 1190 - 1211. 

6. J. Friedel. Metallic alloys. - Nuovo Cimento. 1958. V. 7, 287 - 311. 

7. W. Kohn. Image of the Fermi Surface in the Vibration Spectrum of 
a Metal. Phys. Rev. Lett. 1959. V. 2, 393-394. 

8. W. Kohn and L. J. Sham, Self- consistent equations including ex- 
change and correlation effects. Phys. Rev. 1965, V. 140, A1133- 
A1138. 

9. W. Kohn and L. J. Sham, Quantum Density Oscillations in an In- 
homogeneous Electron Gas. Phys. Review. 1967. V. 137, No. 6A, A 
1697 - A 1705. 

10. W. Kohn and S. H. Vosko. Theory of Nuclear Resonance Intensity 
in Dilute Alloys. Phys. Review. 1960. V. 119, No. 3, 912 - 918. 

11. W. A. Harrison. Solid State Theory. Dover Publ. 1980. 572 pp. 

12. E. M. Lifshits and L. P. Pitaevskii. Physical Kinetics. Oxford: Perg- 
amon (1981). 



16 



13. R. Grassme and P. Bussemer. Friedel oscillations for the tempera- 
tures T^O. Phys. Letters A. 1993. Vol. 175, 441-444. 

14. A. V. Latyshev and A. A. Yushkanov. Longitudinal electric conduc- 
tivity of quantum collisional plasmas / / ArXiv: 1001.3937vl [math- 
ph] 22 Jan 2010. 

15. Mermin N. D. Lindhard Dielectric Functions in the Relaxation-Time 
Approximation. Phys. Rev. B. 1970. V. 1, No. 5, 2362-2363. 



